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SUMMARY 

Aircraft Icing flight research was performed In natural Icing conditions 
with a twin engine commuter type STOL aircraft. In-flight measurements were 
made of the Icing cloud environment, the shape of the Ice accretion on the 
wing, and the corresponding Increase In the wing section drag. Results are 
presented for three Icing encounters. On one flight, the wing section drag 
coefficient Increased 35 percent over the uniced baseline for cruise condi- 
tions while a 43 percent Increase was observed at an aircraft angle of attack 
of 6.2°. 


NOMENCLATURE 

Cj wing section drag coefficient 

C-| wing section lift coefficient 

c wing cord, 1.98 m, 6.5 ft 

KIAS Indicated airspeed, knots 

KTAS true airspeed, knots 

LWC cloud liquid water content, g/m 3 

M Mach number (calculated using aircraft true airspeed) 

MVD cloud median volume diameter droplet size, ym 
OAT static air temperature, °C 

Re Reynolds number (calculated using aircraft true airspeed and wing 
chord length) 

a aircraft angle of attack (referenced to waterline), deg 


INTRODUCTION 

The NASA Lewis Research Center's aircraft Icing flight test program sup- 
ports several program elements. One element of the program Is to develop a 
flight test data base consisting of Icing cloud parameter measurements, Ice 
shapes, and aerodynamic measurements. The data base will be used to validate 



the simulated Icing produced In the NASA Lewis Icing Research Tunnel and to 
validate predictive computer codes. This paper presents In-flight measurements 
of the Icing cloud parameters, the Ice shapes on the wing, and the correspond- 
ing Increases In wing section drag caused by Ice for three natural Icing 
encounters . 

In-flight wing section drag measurements have been performed before for 
airfoils without Ice. (refs. 1 and 2). The experiment consists of measuring 
both total and static pressures across the wake. The pressures are then 
related to the wing section drag coefficient using the Jones equation (ref. 1). 
While this flight experiment has been performed many times for airfoils without 
Ice, no reference was found of It ever being done for airfoils with natural Ice 
accretions on them. The presence of Ice accretions on an airfoil can signifi- 
cantly Increase drag and decrease lift. Many airfoil Icing experiments have 
been performed In wind tunnels (refs. 3 and 4). A major Icing wind tunnel 
facility Is the NASA Lewis Icing Research Tunnel (IRT). The IRT can produce 
simulated Icing conditions through the use of water spray bars and refrigerated 
air flow. Because small differences In Icing parameters can produce large dif- 
ferences In Ice shapes and hence large differences In the aerodynamic effects 
of Ice (primarily at the warmer temperatures (see refs. 4 and 5)) simulation 
techniques should duplicate natural conditions as accurately as possible. 
Therefore a goal of current research Is to compare wind tunnel simulated Icing 
conditions to natural Icing conditions observed In flight. In detail, It Is 
desired to compare (1) the Icing cloud parameters such as cloud liquid water 
content and water droplet sizes (2) the wing Ice shapes resulting from the 
conditions and (3) the change In the wing section drag coefficient. 

To accomplish this task and others the NASA Lewis Research Center Is 
operating an Icing research aircraft and using It to fly Into known Icing con- 
ditions. Work with this aircraft has resulted In three other reports to date 
(refs. 6 to 8). These reports have dealt primarily with Icing Instrumentation 
and overall aircraft performance loss due to Ice. The experiment discussed In 
this report was conducted with the Icing research aircraft during the late 
winter and spring of 1985 and was flown In the Cleveland, Ohio area. 

In addition to the goal of Improving simulated Icing In the IRT, another 
goal Is to use the flight test data to compare with computer predicted results. 
Several computer programs have been developed by NASA to predict droplet tra- 
jectories, Ice accretions, and aerodynamic effects due to Icing (ref. 9). 

While the goals are to make the comparisons described above, this paper pre- 
sents only the flight test measurements of Icing conditions. Ice shapes, and 
drag Increases. 


THE AIRCRAFT 

The Icing research aircraft Is a twin engine commuter type STOL aircraft, 
a DeHavIlland DHC-6 Twin Otter. Maximum gross weight Is 4990 kg (11 000 lb). 
Normal cruise airspeed Is about 259 km/hr (140 KTAS) between the altitudes of 
1829 m (6000 ft) and 3048 m (10 000 ft). Among the various Ice protection 
systems on the aircraft are pneumatic boots on the wings outboard of the engine 
nacelles. The wing section Itself has a NACA 6A series meanllne and a modified 
NACA 0016 thickness distribution (ref. 10). The wing has no sweep back and no 
washout. The wing planform area Is 39 m 2 (420 ft 2 ). 
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INSTRUMENTATION 


Icing Instrumentation 

Measurements were made of cloud liquid water content, median volume drop- 
let diameter size, air temperature, and duration of the Icing encounter. Icing 
environment data was recorded on a digital tape recorder. See figure 1 for 
Icing Instrument locations on the aircraft. 

Instrumentation Included: 

LWC - Johnson and Williams heated wire probe produced by Cloud Technology. 

MVD - Knollenberg Forward Scattering Spectrometer Probe (FSSP) produced by 
Particle Measuring Systems, Inc. (set for a range of 2 to 32 ym) 

OAT - 102CA2W total temperature probe produced by Rosemount Inc. 

It should be noted that the LWC and MVD Icing cloud parameters are diffi- 
cult to measure. Accuracy Is difficult to determine and absolute calibrations 
are lacking. The Instrumentation problem Is very Important, particularly for 
droplet sizes, but It Is beyond the scope of this paper (see refs. 5 and 6). 

An Important fact Is that the above Instruments were tested In the IRT and 
found to be repeatable. Therefore flight versus tunnel comparisons can be 
made. 


Stereophotography System 

A stereo camera system was developed for the Icing research aircraft 
(fig. 2). The system produced stereo pair photo Images of the wing's Iced 
leading edge that enabled measurement of the Ice surface through photogram- 
metric analysis with a minimum acceptable resolution of ±0.0762 cm (±0.03 In.), 
lwo Hasselblad cameras were mounted In the nose of the aircraft behind optical 
glass viewports. The cameras' fields of view encompassed a common portion of 
the wing leading edge located ahead of the wake survey probe. To provide spa- 
tial references In the photo-images needed for photogrammetrlc analysis, an 
array of control points was painted on the wing leading edge section and the 
wing fence. To Increase contrast and character of the Ice surface In the 
Image, a flash unit was mounted on the engine nacelle. The cameras and flash 
were activated from a switch at the copilot's station. The system Is described 
more completely In reference 11. 


Wing Wake Survey Probe 

A wake survey probe was used to measure wing section drag coefficient. It 
was mounted on the wing behind the region where stereophotographs were taken. 
The probe had separate pitot and static probe tips located 1/4 chord aft of the 
wing trailing edge at about 69 percent of the semispan. The probe was driven 
by a motor Inside the wing and could traverse an arc of 180° through the wing 
wake (fig. 3). Static and total pressures were measured by transducers located 
In the cabin. An analog x-y-y plotter In the cabin was used to record the 
pressure traces. During the Icing flights, heaters In the probe kept the probe 
free of Ice and a nitrogen purge system kept pressure lines clear of water. 
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Other Instrumentation 


Altitude and airspeed were obtained from a pitot and static probe mounted 
on a noseboom. The tip of this pitot probe was 2.62 m (8.6 ft) In front of the 
aircraft nose. The position error associated with the noseboom static pressure 
source was calibrated using the trailing cone method. The difference between 
the noseboom and freestream static pressure was approximately 5.9 percent of 
the measured Impact pressure with noseboom pressure greater than freestream. 
This percentage did not vary with angle of attack. The position error correc- 
tion was factored Into the calculation of the aircraft airspeed and altitude. 

An angle of attack sensor was also added to the aircraft. This sensor 
was mounted on the side of the nose on the baggage compartment door. It was 
calibrated during steady level flight against an Inclinometer on a floor rail. 
The angle of attack sensor was made by Specialities Inc. (Model SLZ2303). It 
should be noted that the engineering drawings for the aircraft Indicate that 
the wing Is at a positive 2.5° Incidence to the waterline of the aircraft (the 
waterline Is parallel to the floor line). The angle of attack sensor was ref- 
erenced to the waterline of the aircraft. In addition to the 2.5° Incidence 
of the wing, the local angle of attack of the wing at the section where drag 
measurements were made could differ from that of the aircraft due to the three- 
dimensional flow around the wing. 

On this aircraft, It was unavoidable that the wake survey probe was placed 
behind a region of the wing where ailerons were present. The aileron deflec- 
tions were measured by means of a potentiometer mounted on the aileron control 
cable. 

The altitude, airspeed, temperature, aircraft angle of attack, and aileron 
deflection were recorded manually from digital readouts during each wake survey 
measurement. 


FLIGHT TEST PROCEDURES 

The aircraft was flown In clear air to establish an uniced baseline In 
terms of wing section drag coefficient plotted against aircraft angle of 
attack. This uniced baseline was used as a basis for comparison between the 
Iced versus uniced aircraft. 

The general procedure during Icing flights was (1) accrete Ice while meas- 
uring the Icing environment (2) exit the Icing cloud and document the wing Ice 
shape with stereophotography (3) measure the Increase In wing section drag with 
the wake survey probe at several airspeeds and (4) take additional stereophoto- 
graphs during clear air wake probe traverses to determine If sublimation or 
erosion had altered the Ice shape. 

Airframe Ice was accreted at cruise flight conditions (approximately 
250 km/hr (135 KTAS) ) . While In Icing, a relatively constant cruise airspeed 
was maintained by slowly Increasing power. Icing Instrumentation continuously 
measured Icing cloud parameters during the encounter. The pneumatic deicer 
boots on the wings were not activated during the encounter allowing Ice to 
build up on the wings. After a sufficient amount of Ice had been accreted, 
the aircraft was flown out of the Icing cloud. Generally the quickest means 
was to climb above the Icing cloud. 
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At this point the stereo photographs and wake survey measurements were 
taken. During wake survey runs the airplane was flown steady and level with 
constant airspeed and angle of attack. About three to six measurements were 
usually taken, each at a specified airspeed. Stereo photographs were taken 
for each wake probe measurement. 

During wake survey runs, no aileron Inputs were Intentionally made. The 
ailerons were left In a trim position and the rudder was used If necessary to 
maintain zero degree bank angle. However, the aileron trim position varied 
with airspeed. There was slightly more deflection at the lower airspeeds. 
(See table II). Positive deflections are down. 


DATA REDUCTION 

The Icing Instrumentation data was recorded on digital tape. After the 
flight the tape was read by an IBM 370 computer program which performed neces- 
sary calculations and then plotted the data In corrected engineering units for 
the desired time periods. 

The analysis for the stereo photographs was performed by the Air Force's 
Arnold Engineering and Development Center. From the stereo photographs, the 
x-y-z coordinates of a definitive number of Ice surface points were determined 
along a roughly 0.6 m (2 ft) spanwlse section of the wing. The coordinates 
were then projected onto a single chordwlse plane. Thus the Ice accretion 
profile generated was a "composite" plot (figs. 15, 24, and 33). The scatter 
of the plotted points was Indicative of each accretions roughness and spanwlse 
variability. Reference 11 describes the equipment and data reduction methods 
In more detal 1 . 

The wake survey probe data analysis was done on an IBM 370. First the 
analog plot tracings from the flight tests were digitized and then a program 
was used to reduce the digitized data sets. The program did the following: 
converted raw data to engineering units; applied transducer calibrations; sub- 
tracted reference pressures from the pressure transducer outputs; applied nose 
boom static pressure position error corrections; calculated vertical positions 
of the total and static pressure probes; Interpolated the static pressure data 
In the wake to get static pressures that corresponded to total pressure meas- 
urements; allowed the user to determine limits of Integration; and finally 
Integrated the Jones equation (ref. 2) to determine wing section drag 
coefficient . 

Figures 4 and 5 show typical variations of measured dynamic pressure and 
static pressure across the wake. Note that the static pressure varied consid- 
erably across the wake because the probes were located close to the wing trail 
Ing edge (1/4 chord). The wing has a chord of 1.981 m (6.5 ft) and a max 
thickness which is 16 percent of the chord. The "zero" probe arm position was 
defined as when the probe arm was parallel to the wing trailing edge. In this 
position the probe tips were about 14 cm (5.5 In.) above the wing chord. 
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RESULTS 


Baseline (No Ice) Measurements 

An uniced baseline was established as a basis for comparison between Iced 
versus uniced cases. The baseline consisted of the wing section drag coeffi- 
cient plotted against the aircraft angle of attack. This baseline Is presented 
In figure 6 and represents data from two flights In clear air. The drag coef- 
ficient Is plotted against angle of attack rather than lift coefficient because 
(1) the wing section lift coefficient was not measured during Icing encounters 
and (2) the overall aircraft lift coefficient was altered by the presence of 
Ice accretions on the wing. 

The wake survey probe was also run prior to this season and Initial data 
was published In reference 8. The data acquisition systems, transducers, and 
data reduction programs were entirely different for the two seasons. The 
earlier season employed a computerized data acquisition system that was not 
available for the 1984/85 Icing season. The baseline shown In reference 8 
does not fall on the baseline shown In figure 6. The reason for the discrep- 
ancy Is not completely known, however two reasons justify the publication of 
the new data. The first reason Is that the new data agrees better with numer- 
ical predictions from the Eppler airfoil analysis code (refs. 12 and 13). 

(See fig. 7). The second reason Is that the primary Interest Is In the rela- 
tive effects of the Ice on the drag rather than the absolute numbers. 

Also found In reference 8 Is flight test data showing the uniced wing 
section lift coefficient as a function of aircraft angle of attack. This data 
was acquired by wrapping a pressure belt around the wing and measuring surface 
static pressures. 


Icing Flight Measurements 

Results are reported for three natural Icing encounters for which there 
were adequate Icing conditions and during which research Instrumentation oper- 
ated correctly, providing complete data sets. Data from each Icing encounter 
Is presented separately. For each encounter the order of discussion Is as fol- 
lows: (1) the Icing environment and flight conditions (2) the resulting Ice 

accretion on the wing and (3) the resulting drag Increase caused by the Ice. 
Table I summarizes the Icing cloud conditions for each flight and table II 
summarizes the test conditions for the Individual wake survey measurements. 


Flight 85-17 

The Icing conditions encountered on flight 85-17 caused an accretion 
which was rime Ice with a trough down the stagnation line. The average cloud 
liquid water content was 0.22 g/m 3 . Because of this low LWC the aircraft 
flew In these conditions for 65 min before exiting the cloud. Figure 8 shows 
the large variability of the liquid water content during the encounter. The 
variability Is Interesting considering the fact that most Icing simulations, 
whether In the wind tunnel or on a computer, generally assume constant LWC for 
practical reasons. The average droplet median volume diameter was 12.4 ym. 

The MVD varied little during the encounter compared to the LWC (see fig. 9). 
The average static air temperature was -11.5 °C and was relatively constant 
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(fig. 10). The average true airspeed was 215 km/hr (116 kts) and varied as 
shown In figure 11. The average altitude was 1254 m (4114 ft). Note on 
figure 12 that the aircraft descended Into the Icing cloud at the beginning 
and then descended out of the cloud at the end when enough Ice had been 
accreted. Figure 13 shows the aircraft angle of attack during the encounter. 

Figure 14 presents a photograph of the Ice accretion on the wing. This 
portion of the wing was directly ahead of the area where wake survey measure- 
ments were made. The photographs were originally color and were reprinted In 
black and white. Some of the detail Is therefore lost. Note that the solid 
white color Is characteristic of rime Ice while the trough along the leading 
edge Is a glaze characteristic. 

Figure 15 shows a cross sectional profile of the Ice accretion on the 
leading edge of the wing which was generated from the stereophotography analy- 
sis. As mentioned previously these data points represent locations of the Ice 
surface measured anywhere along a 0.6 m (2 ft) spanwlse section. Thus the 
profile Is a composite Ice shape representing both the general shape and the 
variability of the Ice surface at different locations. The stereo profile Is 
most useful for quantifying the general size of the accretion. 

Figure 16 presents the Increase In wing section drag coefficient due to 
the Ice accretion. Wake Survey probe measurements from the Icing flight are 
compared to the uniced baseline. For the same aircraft angle of attack of 2°, 
the Iced drag coefficient shows an Increase over the baseline of 15 percent. 

A line was fit through these data using the method of least squares. Curve 
fit results are presented In table III. 


Flight 85-24a 

Flight 85-24a was an encounter with mixed Icing (l.e., mix between rime 
and glaze Icing). The average LWC was 0.45 g/m^ (fig. 17). The average MVD 
was 19.5 ym. Note on figure 18 that the FSSP laser probe was not working dur- 
ing part of the Icing encounter. The laser probe Iced up during the first two 
thirds of the encounter. The assumption Is made that the MVD was the same 
during the entire encounter since MVD's have been observed to be relatively 
constant In stratiform clouds. Also the assumption Is made that the laser 
probe was working correctly when It was not Iced up. It Is assumed that It 
was not partially blocked by Ice. From past experience laser probe Icing 
seems to cause total failure of the probe rather than produce merely erroneous 
results. The exception to this Is when large quantities of Ice accrete on the 
laser pod shell and actually change the aerodynamic flow around the probe. 

The amount of Ice on the probe shell can be seen from the aircraft cabin and 
was judged to be acceptable In this case. Considering the assumptions that 
had to be made, the MVD measurement for flight 85-24a was the most uncertain 
of the three flights. 

The average static air temperature was -14.8 °C (fig. 19). The average 
true airspeed was 257 km/hr (139 kts) (fig. 20). The average pressure alti- 
tude was 2008 m (6588 ft) (fig. 21). The average aircraft angle of attack was 
0.9° (fig. 22). 

Figure 23 shows a photograph of the Ice accretion on the wing. The Ice 
Is nearly clear In the center which Is a glaze characteristic. Along the edges 
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of the accretion the Ice Is more like rime. One might wonder why flight 85-17 
was rime Ice and flight 85-24a was mixed Ice when the temperature was colder 
on flight 85-24a. The reasons are that flight 85-17 had a much lower LWC and 
smaller droplet sizes so the freezing fraction was greater. Figure 24 pre- 
sents the profile generated from the stereo photographs. 

Figure 25 shows the Increase In wing section drag coefficient due to the 
Ice accretion. At an aircraft angle of attack of 2°, the drag coefficient 
Increased 23 percent over the baseline. Six measurements were taken with Ice 
on the wing. It Is believed that measurements 4 to 6 are lower than they 
should be due to Ice erosion and sublimation; therefore they were not Included 
In the figure. The change In the Ice shape was caused by the fact that after 
sufficient Ice had accreted on the aircraft It was flown out of the clouds, 
Into the sun, to terminate Icing. The sunlight hastened the sublimation of 
the Ice even though the temperature was still cold. The general procedure was 
to fly away from the sun, (when possible) In an attempt to shield the Ice 
accretion from the sun for as long as possible. Also to minimize sublimation, 
measurements were begun with the low speed data points first and then followed 
with the higher speed data points. 


Flight 85-24b 

Flight 85-24b was another mixed Icing encounter. The aircraft was In 
Icing conditions for 20 min. The average LWC was 0.46 g/m 3 (fig. 26). The 
average MVD was 15.1 vm (fig. 27). Note that the laser probe Iced over 
towards the end of the encounter when the LWC and Icing rate were very high. 

The average static air temperature was -14.3 °C (fig. 28). The average true 
airspeed was 259 km/hr (140 kts) (fig. 29). The variations In airspeed were 
caused by turbulence In the Icing cloud. Also, the pilot sometimes changed 
altitude slightly to stay In the best Icing conditions. The average altitude 
was 1976 m (6483 ft) (fig. 30). The average aircraft angle of attack was 1.2° 
(fig. 31). 

A photograph of the Ice accretion on the wing Is shown In figure 32. The 
corresponding profile generated from stereophotography analysis Is shown In 
figure 33. The accretion Is very similar to that of Flight 85-24a since the 
environment was nearly the same, but this accretion Is slightly larger; with 
small "horns" that were probably caused by the additional five minutes In Icing 
conditions . 

The resulting Increase In wing section drag coefficient over the uniced 
baseline Is presented In figure 34. Due to the larger Ice accretion, the drag 
Increase Is greater for this case. At a 2° angle of attack, the drag Increase 
was 39 percent of the uniced baseline. 


CONCLUSIONS 

Three data sets have been collected which can be used to begin compari- 
sons between natural Icing In flight and simulated Icing; either In an Icing 
research tunnel or using computer codes. Flight 85-24b Is the best data set 
for comparison purposes because the largest Ice accretion and the largest drag 
Increase were observed on that flight. Flight 85-24a was very similar to 
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85-24b except the Ice accretion was slightly smaller and the drag rise less. 
Flight 85-24a had the most uncertain MVD measurement of the three encounters. 

Flight 85-17 had the most data points associated with the drag plot but 
unfortunately the drag rise was rather small. Also the LWC was very low 
(0.22 g/m^), and the Icing duration was very long (65 min) making the 
encounter difficult to simulate. 

All three flights had highly variable LWC's during the encounter. This 
variability could be compared to the constant Icing conditions which are pro- 
duced In an Icing wind tunnel. All three flights had similar Ice shapes which 
are a mix between rime and glaze Ice. This Is primarily due to the tempera- 
tures being relatively close. It would be highly desirable to document addi- 
tional encounters which produce more pure rime and glaze Ice accretions (l.e., 
Ice formed at very cold temperatures and at temperatures very near freezing, 
respectively). Also, flights are desired where the accretions are large. A 
larger Ice accretion takes on a more unique shape which Is a function of the 
physics Involved. These larger shapes are then a more significant test when 
comparing flight results to simulated Icing results, either In the wind tunnel 
or on the computer. 
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TABLE I. - SUMMARY OF ICING CLOUD CONDITIONS 


Flight data 

F1 1 

ght number 


85-17 




Date 



EHBSHOSlil 

4/2/85 

1. Average pressure altitude, m 

1254 

2008 

mm 

2. Average true airspeed, km/hr 

215 

257 

mm 1 

3. Average aircraft angle of 
attack, deg 

2.4 

0.9 

1 .2 

4. Average right aileron 
deflection, deg 

-1.2 

-1.1 

-0.9 

5. Start time of Icing encounter 

09:48 

09:55 

10:49 


Icing cloud data 


6. Average static temperature, °C -11.5 -14.8 -14.3 

7. Average cloud liquid water 

content, g/m^ 0.22 0.45 0.46 

8. Average median volume droplet 

diameter, ym 12.4 19.5 15.1 

9. Duration of encounter, min 65 15 20 

10. Extent of encounter, km 233 64 86 


Ice accretion properties 

1 1 . Type of Ice 

12. Shape of Ice (see photo and 

profile) 

Rime 

Mixed 

Mixed 

















TABLE II. - TEST CONDITIONS AND WING SECTION DRAG DATA 














TABLE III. - CURVE FITS OF DRAG DATA 
(LEAST SQUARES METHOO) 
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Figure 5. - Typical variation of static pressure across wing 
wake * flight 85-12, run 9. 
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Figure 6. -Variation in wing section drag coefficient with aircraft 
angle of attack for wing with no ice. 



Figure 7. - Comparison of flight test baseline results with 
numerical prediction (for wing with no ice). (Prediction 
assumes that winq angle of attack equals aircraft angle 
of attack plus 2. 5 . ) 
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Figure 9. - Cloud median volume droplet diameter during icing 
encounter - flight 85-17. 
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Figure 10. - Static air temperature during icing encounter 
flight 85-17. 
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Figure 11. - True airspeed during icing encounter - flight 85-17. 
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Figure 12. - Pressure altitude during icing encounter - flight 85-17. 



Figure 13. - Aircraft angle of attack (referenced to aircraft waterline) 
during icing encounter - flight 85-17. 



Figure 14. - Ice accretion on wing - flight 85-17. 
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Figure 19. - Static air temperature during icing encounter - 
flight 85~24a. 



Figure 20. - Aircraft true airspeed during icing encounter - flight 85-24a 
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Figure 24. - Ice accretion profile generated from stereo 
photography - flight 85-24a. 
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Figure 25. - Increase in wing section drag coefficient due to mixed 
ice - flight 85-2 4a. 
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Figure 26. - Variation in cloud liquid water content during icing 
encounter -flight 85-246. 
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Figure 27. - Cloud median volume droplet diameter during icing 
encounter - flight 85-24b. 
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Figure 28. - Static air temperature during icing encounter - 
flight 85-246. 
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Figure 29. - Aircraft true airspeed during icing encounter - flight 85-24b. 
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